The Ferromagnetic Heisenberg XXZ chain in a pinning field 
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We investigate the effect of a magnetic field supported at a single lattice site on the low-energy 
spectrum of the ferromagnetic Heisenberg XXZ chain. Such fields, caused by impurities, can modify 
the low-energy spectrum significantly by pinning certain excitations, such as kink and droplet states. 
We distinguish between different boundary conditions (or sectors), the direction and also the strength 
of the magnetic field. E.g., with a magnetic field in the z-direction applied at the origin and ++ 
boundary conditions, there is a critical field strength B c (which depends on the anisotropy of the 
Hamiltonian and the spin value) with the following properties: for B < B c there is a unique ground 
state with a gap, at the critical value, B c , there are infinitely many (droplet) ground states with 
gapless excitations, and for B > B c there is again a unique ground state but now belonging to 
the continuous spectrum. In contrast, any magnetic field with a non-vanishing component in the 
xy-plane yields a unique ground state, which, depending on the boundary conditions, is either an 
(anti)kink, or an (anti)droplet state. For such fields, i.e., not aligned with the 2-axis, excitations 
always have a gap and we obtain a rigorous lower bound for that gap. 

PACS numbers: 75.10.-b, 75.60.-d, 05.50.+q 
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I. INTRODUCTION 

The quantum spin-j Heisenberg XXZ chain has been 
the focus of intensive studies in recent years. The spin 
1/2 chain by itself has connections with a surprising vari- 
ety of interesting mathematical structures, such as quan- 
tum groups 1 , vertex algebras 2 , and fundamental prob- 
lems in combinatorics 3,4 , to name just a few. Of course, 
the interest in the XXZ model is not limited to mathemat- 
ics. In 1995, Alcaraz, Salinas, and Wreszinski 5 and, inde- 
pendently, Gottstein and Werner 6 discovered that, with 
suitable boundary terms, the ferromagnetic XXZ chain 
possesses a family of kink ground states which describe 
a domain wall of finite thickness (these domain walls are 
exponentially localized, with a width depending on the 
anisotropy parameter A, which diverges as A J, 1). More- 
over, it was shown in 5 , that similar states exist for the 
spin j model for arbitrary j and in all dimensions. 

From the physical point of view, the discovery of Gi- 
ant Magnetoresistance and its connection with transport 
properties in the presense of magnetic domain walls has 
also spurred renewed interest in the microscopic descrip- 
tion of domain walls 7-16 . Of particular relevance are 
low-lying excitations associated with them. Koma and 
Nachtergaele discovered that, although the XXZ model 
has a gap in its spectrum above the trivial translation 
invariant ground states, gapless excitations exists associ- 
ated with diagonal domain walls (11, 111, ...) in two or 
more dimensions 17 . The scaling behavior of these exci- 
tations was recently determined in 18 and , . 

It is interesting to note that the kink and antikink 



ground states were discovered by a careful study of the 
XXZ chain with the special boundary conditions that 
make the spin 1/2 model S , C/ g (2)-symmetric. Although 
this quantum group symmetry is destroyed for j > 1/2 
or d > 1, interface ground states exist in general. In one 
dimension it has been proven rigorously that no other 
ground states exist, in the sense of local stability, no mat- 
ter what boundary conditions are considered 21,22 . In the 
last reference it was also proved that the XXX chain does 
not have domain wall ground states that are stable in the 
infinite- volume limit. 



There is an obvious need for a clear understanding of 
excitations near magnetic interfaces in order to develop 
more accurate models of electron scattering at such in- 
terfaces. It has been noted, however, that pinning of 
interfaces by impurities may have to be taken into ac- 
count as well 9 . Here, we study the XXZ chain perturbed 
at one site by a magnetic field as a caricature model for 
a pinned domain wall. Admittedly, the one-dimensional 
nature of the model restricts its direct applicability to 
experimental situations. We will see however that the 
low-lying spectrum of this model of pinned interfaces ex- 
hibits a number of interesting features that we expect 
will carry over, mutatis mutandis, to the two- and three- 
dimensional case. 



Let us now define the model precisely and briefly sum- 
marize our main results. The spin j XXZ Hamiltonian 
without boundary terms and with anisotropy parameter 
A > 1 is defined on the finite chain labeled by the inte- 
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gers from a to b by 
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where the copies of the spin operators at position x, 
S%, S% satisfy the usual commutation relations: 



f qa c/3i _ ■ a/37 C7 
l°x J °x\ — 1 t °x- 



(2) 



The anisotropy A > 1 has been put in front of the XX 
part so that we can easily take the Ising limit A — > oo. 
Unless otherwise stated, we set a = 1. We consider the 
perturbation of Hq obtained by adding a term B ■ S y , 
i.e. a magnetic field at the site y. As a way to impose 
boundary conditions, we also add magnetic fields in the 
z-direction at the boundary spins. First, consider bound- 
ary fields in the negative z-direction at both ends, which 
we will refer to as (H — h) boundary conditions (b.c), in- 
dicating that they favor the spins at the ends to point in 
the positive z-direction. 

If the perturbation at the interior site y has a com- 
ponent orthogonal to the z-direction, we find that the 
ground state is unique and describes a droplet state, i.e., 
the magnetization is reduced from its maximum possible 
value in some neighborhood of y. Strictly speaking, the 
magnetization is reduced everywhere in the chain, but by 
an amount that decays exponentially fast away from y. 

However, when the field is in the z-direction, then there 
is a critical value B c such that for B < B c , the all spin-up 
state is the ground state. At B = B c , there are infinitely 
many ground states which are droplet states describing 
domains of negative magnetization of arbitrary size em- 
bedded in a environment of positive magnetization. For 
stronger fields, B > B c , the magnetic field selects the all 
spin-down state as its ground state. This is illustrated in 
Figure 1. 
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FIG. 1: Ground state energy (thick line) in a magnetic field 

in the z-direction and (H — h) or ( ) b.c. The dotted lines 

indicate droplet states which describe excitations except at 
one value of the field strength, B c , given in Proposition II. 4. 
Also see Figure 6. 

The ground state picture is simpler when we impose 
(H — ) b.c, i.e. fields in opposite directions at the bound- 



ary spins. For boundary fields of magnitude (3), and 
without a perturbation in the interior, we then have a set 
of kink states as the ground state, one for every possible 
value of the magnetization 5 ' 6 . In that case, any non-zero 
field B at an interior site y selects a unique ground state. 
If B is parallel to the z-direction the ground state is in 
the continuous spectrum, and hence there are excitations 
of arbitrary small energy. If there is a non- vanishing com- 
ponent of B in the xy-plane, the unique ground state is 
separated by a gap form the rest of the spectrum. The 
unique ground state is a kink state centered at a position 
which we calculate. We also obtain an estimate for the 
gap. 

In Section 2, we define the model and find the set of 
ground states. Section 3 is devoted to the study of the 
gap in the spin 1/2 case. Some less illuminating calcula- 
tions are presented in three appendices. 



II. THE MODEL AND ITS GROUND STATES 

The main lesson to be learned from the proof of com- 
pleteness of the list of ground states of the infinite fer- 
romagnetic XXZ chain (cf. 23 ), is that one only needs 
to study finite chain Hamiltonians with very simple b.c. 
This remains true if we add a bounded perturbation with 
finite support to the Hamiltonian, e.g., a magnetic field 
at one site. These simple b.c. are fields in the z-direction 
with either equal or opposite sign which we will introduce 
shortly. 

The Hamiltonian defined in (1), is non-negative, 
and the two translation invariant all spin- up/down states 
are the ground states of Hq. It will be convenient to 
separate them by adding the equal-sign boundary fields, 
jA(S 3 + S 3 ) with 



(3) 



A = A(A) = \f\ - A" 2 , 

and define the droplet and antidroplet Hamiltonian 

= H -jA(S 3 + S 3 -2jl), (4) 
-jA(S 3 + S 3 + 2jl). 



H - H Q 



For convenience we have normalized the ground state en- 
ergy to 0. By reflecting all S 3 into— S 3 the two Hamil- 
tonians, H ++ and H are unitarily equivalent and we 
only study Hq + . 

Additional ground states emerge when we add 
opposite-sign boundary terms. It turns out that precisely 
for the fields ±jA(S 3 — S 3 ) one discovers the full set of 
new ground states. Therefore, we define the kink and 
antikink Hamiltonian 



H- 



H -jA(S 3 -S 3 ), 
H +jA(Sf-Sl). 



(5) 
(G) 



Again, by spin reflection, the kink and anti-kink Hamil- 
tonians are unitarily equivalent. Let us define 



~k(^x^x+l + Sx^x+l) Sx^x+l +J 2 1 
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-jA(S x S 3 X+1 ), (7) 

and 

+.M(S 3 -S 3 +1 ). (8) 
In terms of these interactions terms we may write 

x=l 

and 

6-1 

ff " + = Etti- ( 10 ) 

x=l 

This will be used in Section 3. 

It is useful to introduce another parameter, < q < 1, 
such that q + q^ 1 = 2 A. The Hamiltonians defined in 
(5) and (6) in the spin 1/2 case, commute with a repre- 
sentation of SU q (2). For j > 1/2, the only obvious con- 
served quantity is the total S 3 -component, which com- 
mutes with all the Hamiltonians defined above. 

In the following we first deal with the kink Hamilto- 
nian, Hq~ . There is a unique ground state for each value 
(sector) m of S 3 , all of which have the same energy 0. The 
eigenvalues, m x , of S 3 are in {-j, -j + 1, . . . , j - 1, j}, 
such that the total S^-component takes the values m = 
J2xe[i 6] m x- The eigenvectors of S 3 are denoted by 
\(m x )}, and we have Sy\(m x )) = m y \(m x )}. Further, 
let 

The unique ground states in the respective 5* 3 -sectors are 
called kink states which were found by Alcaraz, Salinas 
and Wreszinski 5 . They are given by 

b 

^=EI1 q- xU - m ^w m J(m x )), (11) 

{m x } x=l 

As the ip m are not normalized, we also define 4> m = 
V'm/||V ; m|| ■ The sum over {m x } is restricted to combi- 
nations such that "^2 m x = m. If m = ±jb is maxi- 
mal/minimal, then the state 4> m is the all spin- up/down 
state, i.e. the magnetization profile in the z-direction, 
{4>m\S x \4> m ) = ±j- Both from a physical and mathemat- 
ical point of view, the infinite chain limit is the most 
interesting case. Clearly, more care has to be taken 
when using an infinite volume Hamiltonian. The nat- 
ural tool is the GNS representation. In our case here, 
there are four (unitarily) incquivalcnt representations on 
Hilbert spaces, which are also called sectors. These 
are the (anti)kink sectors containing the infinite volume 
(anti)kink ground states, and the (anti)droplet sectors 



with the all spin (down)up ground state. As mentioned 
above, for the infinite volume limit it is sufficient to study 
the effect of the perturbation for the finite chain Hamil- 
tonians, Hq + and Hq~. For more details, we refer to 23 . 

For our purposes it will be very convenient to define 
the states 

1/>(Z)=C E ^^^m. (12) 
\m\<jb 

where C is a normalization constant. They are product 
states, i.e. 

6 

#*) = ®X*(*), (13) 

x=l 

with 

Xx (z) = (1 + \z\ 2 q- 2x )-3 E (zQ- X y- m *™m x \(m x )). 

(14) 

The same construction can be carried through for the 
antikink Hamiltonian. We denote the corresponding 
states for the antikink Hamiltonian by tp(z). If we let 

**(*) (15) 

= (l + \z\ 2 q- 2x yi ]T (zq- x ) j - m ^mj(-m x )), 

m x =-j 

then 'ijj(z) = ® x=1 Xx{ z ) are ground states of H h . 

Now, let B — (Bi, B2, B3) be a magnetic field vector 
with (real) parameters, and V = B ■ S = BxS 1 + B 2 S 2 + 
B 3 S 3 . Then the eigenvalues of V are ||S|| • m with m — 
-3,-3 + !. ■ • ■ ,3- Define 

H+-(B) = H+-+B-S y , (16) 
H++{B) = H++ + B-S y . (17) 

In the study of the spectrum of H A (B) and H ^(B), 
it is important to distinguish two cases: B\ + B\ > 
and B\ + B\ = 0. The ground state in these two cases is 
described in Propositions II. 1 and II. 3, and the Remarks 
II.2. 

PROPOSITION II. 1 (Kink sector, B\ + B\ > 0). 

Let 1 < y < b. Then, the ground state of H + ~(B) is the 
state ip{z) of eq (13) with z = — |f l~jg 3 q y ■ Its energy is 
-3\\B\\. 

The proof is a combination of previously known results 
and a straightforward calculation. Sec Appendix A. 

REMARKS 11.2. 0. The ground state of H~+(B) 
is the state ip(z) of equation (15) with 5 = 

— ^gj'+igl q v ■ This follows by a rotation by the an- 
gle tt with respect to the x-axis. 



1. For simplicity, let us assume that B\ + B\ = 1. 
If B 3 = 0, then the ground state, tp(z — — + 
1B2) q y ), is a kink state (exponentially) localized at 
the magnetic field at y; among the spanning set of 
ground states ^{z), the perturbation picks the one 
which is most localized at y. If B 3 7^ 0, the extra 
term B 3 Sy has the effect of shifting the kink from 

y by the distance | \og q (i/l + B 3 + B 3 )\ to the left 
if B 3 > 0, and to the right (by the same distance) 
if B 3 < 0. 

2. The proof also shows that the state i^(z) with 
z = + ^Lt_Ba 1 V ^ s an e ig ensta te of H + ~(B) with 
energy +j||_B||. We see this (it is obviously linear if 
the field is in the x or y direction only, see Figure 
2) branch ascending from I-B3I/2 in Figure 2-3. In 
Figure 4 it has too high an energy to be among the 
plotted lowest eigenvalues. 

3. Of special interest is the second-lowest eigenvalue, 
in particular, whether there is a gap uniformly in 
the number of sites, 6, and how it depends on B and 
the anisotropy A. This will be treated in Section 3. 
We can extend a method 2 which was first applied 
to prove a gap for Hq . 

4. We discuss now qualitatively the low energy spec- 
trum, and assume for simplicity that j = 1/2 and 
B 2 = 0. It was proven in 23 that the gap above the 
6+1 ground states of the unperturbed Hamiltonian, 

, is equal to 1 — cos (7r/6)A _1 , which tends to 
1 — A -1 in the infinite chain limit. 

If B 3 = 0, then there are b eigenvalues (re- 
call, b is the length of the chain) of (B,0, 0) 
descending from 0. This can be seen as fol- 
lows. We introduce the function N(B) = 
X(-oo.o] (H + ~ (B,0, 0) — \B\) counting the num- 
ber of non-positive eigenvalues; here X(-oo,o] is 
the characteristic function of (— 00, 0]. N(B) is 
monotonically increasing, and equal to b + 1 for 
\B\ < 1 - A" 1 + 0(fr _1 ); this is guaranteed by 
the gap above the ground states. In item 6 below 
we calculate the lowest energy state, ip e , descend- 
ing from 1 — A -1 at B = with energy equal to 
1 — A" 1 - ||_B| (uptoO(&- 1 )). This state is 'paral- 
lel' to the ground state energy and intersects with 
the state in item 2 at \B\ = 1 — cos (7r/6) A -1 , see 
Figure 2. 

We can say more about the average of these lowest 
6+1 eigenvalues by recalling the Min-max Princi- 
ple, namely that their average is a concave func- 
tion in B. By symmetry (B — > — B) its (left and 
right) derivative is always negative and less than 
1/2. Since the average is at B = it continues 
to be negative. In the infinite volume limit there 
is an infinite number of ground states of , and 
eigenvalues for H A (B, 0, 0) have to accumulate at 
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FIG. 2: Energy of the lowest 16 eigenvalues for the spin 1/2 
kink Hamiltonian for the field (£?, 0, 0) on a chain of 13 sites 
and A = 2.25. Notice the ground state energy (straight line 
downwards from 0), the gap to the 2nd eigenvalue, the energy 
and the energy 1 - A -1 - \B\/2 of ip e . 

some value between — \\B\ and 1 — A -1 — \\B\, see 
Figure 2, and in case B 3 ^ 0, see Figures 3-4. 

5. If B 3 + 1 0, then there appear to be b — y + 1 eigen- 
values descending from —\B 3 \/2; b — y + 1 is the 
number of kinks to the left of y, see Figure 3. The 
next eigenvalues depend on B 3 . H \B 3 \ < I — A -1 , 
then the next lowest y eigenvalues descend from 
-\B 3 \/2. If \B 3 \ > 1 - A -1 , then the state ip e ap- 
pears to be the next lowest in energy, see Figure 
4. 

6. We can calculate the state mentioned in the previ- 
ous two items. Let ifj e = J2x=i a ^^IT) be the first 
excited state of Hq in the one-overturned spin- 
sector. The coefficients a x are a solution to the 
discrete Laplace equation (see 22 , or cf . Appendix C 
by setting B = 0). The energy of ip e is equal to the 
gap 1 - A" 1 + Oib- 1 ). Now, define 

b b 

71— X — l 

By choosing z as in Proposition II. 1, we obtain the 
equation 

H+-{B)^ e {z) = (l - A" 1 - ±11511) 4> e (z) + O (6- 1 ) . 

In Figure 2, this is the straight line parallel to the 
ground state energy. 

Next, we consider the kink sector and magnetic fields of 
the form B x = B 2 = 0,B = B 3 ^ 0. Let B > 0. Then, as 



-1.5 



FIG. 3: Energy of the lowest 16 eigenvalues for the spin 1/2 
kink Hamiltonian for the field (B, 0, A/6) on a chain of 13 sites 
and A = 2.25; the ^-component, B3 — A/6 — \/l — A _1 /6 is 
chosen small compared to the gap 1 — A -1 of Hq . Notice 
the ground state energy, the gap above it, and the branches 
descending from ±£3/2, and 1 — A -1 — B3/2. 



b, the size of the system increases, the ground state tends 
to the all spin-down state, 1 1). This vector is no longer in 
the infinite volume kink-sector very much as e lkx is not a 
genuine (i.e., normalizable) eigenvector of the Laplacian 
on the real line. In other words, | J.) is part of the continu- 
ous spectrum. So let us consider the orthogonal sequence 
of kink states, ip n ; n, as usual, is the total z-component. 
Then, the sequence {ip ni H^ (0, 0, B), tp n ) converges to 
—jB as n — > —00. Since the spectrum is closed and —jB 
is the least possible eigenvalue it has to be the ground 
state energy. Therefore, in the infinite chain limit, — jB 
is contained in the continuous spectrum, and is hence 
non-isolated. We conjecture that there is no other contin- 
uous spectrum close to — j\B\, and thus — j\B\ is purely 
an accumulation point of eigenvectors. We do not give 
a proof of this here. Similarily, if B < 0, then the bot- 
tom of the spectrum is jB and there is no gap above the 
ground state, which is obviously the all spin-up state. We 
illustrate the low energy spectrum in Figure 5. 

Let us collect our results in the following proposition: 

PROPOSITION II.3 (Kink sector, B\ + B\ = 0). 

The bottom of the spectrum of (0,0, B) is equal to 
—j\B\, which is part of the continuous spectrum. Excita- 
tions above the ground state are gapless. 

Now we consider the Hamiltonian H ++ (B). It is use- 
ful to decompose this as a sum of a kink and anti-kink 
Hamiltonian (b > 3): 



FIG. 4: Energy of the lowest 20 eigenvalues for the spin 1/2 
kink Hamiltonian with the field (B, 0, 3) on a chain of 13 sites 
and A = 2.25; the field in the z-direction is chosen large 
compared to 1 — A -1 . Notice that the branches stemming 
from B3/2 are too high in energy to be plotted here. We see 
the ground state energy, the gap above it, and the energy of 
ipe descending from 1 — A -1 — -B3/2. In between there are 
b — y + 1 = 6 states bending downwards from — -B3/2 = —1.5. 

and thus 

H++(B)-2j 2 Al = H+-(f ,f ,f -jA) 

+ H M& 'f >f -MM") 

We start with the case Bf + B\ > 0. As in the kink 
sector we will find a unique ground state. Let 

\\(B 1 ,B 2 ,B 3 -2jA)\\+B 3 -2jA 

z = a" , 

B 1 -iB 2 q ' 

then according to Proposition II. 1, y i(,z) is the unique 
ground state of (-Bi/2, B 2 /2, B 2 /2 - jA), while the 

anti-kink state ip\ y u(z) is the corresponding ground state 
of Hr+(Bx/2, B 2 /2, B 3 /2 - jA). They are both prod- 
uct states which happen to satisfy Xy( z ) = (~ l) 2j Xy C&)> 
because 

\{B x Sl + B 2 S 2 y + (B 3 - 2j A)S 3 y )x i (z i ) 
= -l\\(B 1 ,B 2 ,B 3 -2 J A)\\x i (zt), 

where x ( z ) stands for either Xy( z ) or Xy{ z )- Thus 
y b 

x—1 x— y+1 

is the unique ground state of Hq + (B) with energy 
— ill (Bi,B 2 , B 3 — 2jA)\\ + 2j 2 A. 
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FIG. 5: Low energy spectrum in the spin 1/2 kink sector 
for the field (0,0, 1.5) on a chain of 11 sites with A = 2.25. 
The high energy spectrum is similar to Fig. 7 and will not be 
reprinted, n is the total z-component. 



FIG. 6: The five lowest eigenvalues in the spin 1/2 droplet 
sector for a chain of 13 sites with magnetic field (0, 0, B) and 
A = 2.25. The eigenvalues are numbered by n, the number 
of overturned spins; e.g. n — means the all spin-up state. 



Similar to the kink-sector, the vector 
v b 

x—1 x— y+l 

is another eigenstate with energy j\\ (Bi, B2, B^ — 2jA)\\ 4- 
If A. 

Finally, we come to the case Bx = B% = 0,B = B 3 . 
When B = 0, it was proven in 25 that for spin 1/2, and 
in the infinite volume limit 1 — A -1 is the gap above the 
ground state. It can also be shown that there exists a 
gap, S, for higher spins although no precise estimates are 
known. This implies that uniformly in the size of the 
lattice, that the all spin-up vector is the unique ground 
state for B < B c = <5/(2j), where S is the (strictly posi- 
tive) gap of Hq + . As we mentioned in the Introduction, 
the value B = A/(2j) is very particular and interesting. 
It has been analyzed 26 in the context of droplet states 
for spin 1/2 but again the method extends to general j. 
In fact, the set of ground states is infinitely degenerate 
(in the infinite volume) and consists of pairs of symmet- 
ric kink-antikink states (i.e. droplets), all of which have 
the same energy jA. The magnetization profile in the 
z-direction is symmetric with respect to the center of the 
field at y. Excitations are gapless because large droplets 
which are antisymmetric with respect to y come arbitrar- 
ily close in energy to jA. 

Since the ground state energy is concave, and since 
for B = and B — A/2j, the all spin-up vector is 
a ground state, we conclude that for all B < A/(2j), 
the all spin-up vector is the unique ground state. Nu- 
merical experiments for spin 1/2 indicate that in the 
region B < A the eigenvalues are ordered by their to- 



tal S^-value such that the second-lowest eigenstate is 
in the sector with one overturned spin, and has energy 
8 (B) = 1 - VA~ 2 +B 2 + see (C3); its (infinite 
volume) derivation is given in Appendix C. The third- 
lowest eigenvector is in the sector with two overturned 
spins, and so on. The lowest eigenvalues with respect to 
the total S^-component accumulate at the line A — B/2. 
They all meet at the critical value A, where the ground 
state becomes infinitely degenerate. Assuming that this 
ordering holds true, we conjecture that the gap for spin 
1/2 and B < A equals 1 - \/A- 2 + B 2 + - B), 

which converges to 1 — A for large —B, and vanishes at 
B = A. 

For B > A/(2j), the all spin-down state is the unique 
ground state with energy A/(2j) — jB, which is part of 
the continuous spectrum; in fact, A/(2j)—jB is purely an 
accumulation point of eigenvectors. It seems that for B > 
A/2j the eigenvalues are also ordered according to their 
total z-component, n, but this time in the opposite way. 
I.e. lower n means lower energy, and clearly the lowest 
is the all spin-down state. Similar to the kink sector, 
we will not prove here that the rest of the continuous 
spectrum is separated from the ground state. 

The situation is illustrated in Fig. 6-8 for j = 1/2, 
cf. also Figure 1. Let us summarize our results in the 
following proposition: 

PROPOSITION II.4 (Droplet sector). The ground 
state of the droplet Hamiltonian, H ++ (B), on a chain 
of length b > 3 depends on the magnetic field B in the 
following way: 

1. If B\+B\ > 0, then the ground state is unique. The 
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Spectrum of H ++ (0,0,B) for L = 13, B=1.5A 



these ground states are gapless. For B > B c , the all 
spin-down state is the unique ground state, which 
is an accumulation point of eigenvectors. Hence, 
excitations are gapless. 



■ 

I 



10 11 12 13 14 



FIG. 7: Here, we plot the full spectrum of the spin 1/2 droplet 
Hamiltonian H ++ (0,0, 1.5A) for 13 sites, and A = 2.25. The 
index n is the number of overturned spins. The ground state 
is the all spin-down state, i.e. n = 13. 



Amplified low energy spectrum of H ++ {0,0,B) 



FIG. 8: Here, we amplify the low energy spectrum for the 
same Hamiltonian as in Figure 7. At n — 13 we have indicated 
the ground state. Clearly, one sees monotonicity of energy vs 



ground state energy is 
2fA. 



-j\\(B 1 ,B 2 ,B 3 -2jA)\\ 



2. If Bi = B 2 = 0,B = B 3 , and B c = A/(2j), then 
for B < B c the unique ground state of H ++ (0, 0, B) 
is the all spin-up vector with energy jB. For B = 
B c the ground states are droplet states which are ( in 
the thermodynamic limit) infinitely degenerate with 
energy jA. In infinite volume, excitations above 



III. ESTIMATE FOR THE SPECTRAL GAP IN 
THE CASE 3 = 1/2 

Here we prove a uniform lower bound on the difference 
between the ground state energy and the energy of the 
first excited state for the spin 1/2 Hamiltonians H^AB) 

and H^~t(B) on a finite chain [1,6] with the impurity 
field at y. 

Before we prove these gap inequalities we introduce 
the methods which were invented in 24 and 27 . Let 
Cj, i = 0, . . . , N be a sequence of connected intervals with 

UiI=o ^» = [1' ano - such that two intervals have at most 
one lattice point in common. Let hc t > be some (local) 
Hamiltonians acting on C 2 ' Ci ', and define 



H [1M 



N 

E 

i=0 



(19) 



H[x,b] ac ts on Hb = ® i=1 C 2 . We assume that 
kerif[i.b] 7^ {0}. Let 74 denote the gap of hc i: i.e. the 
smallest non-zero eigenvalue of /ic 4 ■ It is clear that 



ker#[i >6 ] = P| ker he- 



rn 



i=0 



Let A C [1,6], then we define G\ to be the orthogonal 
projection onto the 



ker 



E 

:CiCA 



ha 



(21) 



We use the convention that if Cj (f_ A for any i, then 
we set Ga = 1. From these definitions we derive the 
following properties: 

1. G A G A , = Ga-Ga = G A < if A c A'. 

2. G A G A , = Ga-Ga if A n A' = 0. 

3. h Ci >n{i-Ga)- 

Next we define the intervals Aj = Uj<; Pjj an d operators 
Ei,i> 0, on n h by 




for i = 

for 1 < i < N 

for i = N 



(22) 



These operators are mutually commuting projections 
adding up to 1, i.e. 



A 



E* = Ek, EiEj = SijEi, E i = 1 - 



(23) 



i=0 



The key assumption in order to deduce a gap for H\im 
from the gaps of ha is the following assumption: 
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ASSUMPTION III.l. There exists a positive constant e 
such that < e < 1/ y/2 and 

E l G Cz+1 E l < e 2 E t , < i < N - 1, (24) 

or equivalently, 

\\G Ci+1 Ei\\ < e, < i < JV- 1. (25) 

Further, we assume that the gaps, 7^, for the local Hamil- 
tonians are bounded from below, i.e. 7i > 7 > 0. 

The conditions (24) and (25) are equivalent due to 
Gc z+l Ei — Gc z+l G\ z — Ga 4+1 . 

Now we are ready to state the main theorem which we 
apply in all three case below. 

THEOREM III.2 (Nachtergaele 24 ). With the above 
definitions and under the assumptions in (25) let ij) be 
orthogonal to the ground states of • Then 



(V>,i? [M] V)>7(l-V2e)< 



(26) 



I.e. the gap in the spectrum of H[i.h\ above is at least 
7(1 -V2e) 2 . 

Proof. Let tp be orthogonal to the ground state, i.e. 
G [ltb] il> = 0. Then, ||V>|| 2 = £ <„ <JV \\E n yj\\ 2 .^ 

We estimate ||-E„-0|| 2 in terms of (f/>, Hc n (B)tp) as fol- 
lows. First notice, that for m < n — 2, or m > n + 1, 
E m Gc n = Gc n E m . Now we insert Ge„ and the resolution 
of identity, {E n } 7 and we get 

\\E n n 2 

= (V>,(l-G c J£ n V) + (V>, J2 E ™G Cn E n ^) 

0<m<N 

= (V>, (1 - G c jE n ^) + ((£„_! + E n )$, G Cn E n ^). (27) 
Let ci, c 2 > 0, then 



\E n i>\\ 2 

< -L^,(l-G Cn )^) + ^,E n i>) 



(28) 



1 



C 2 , 



+ — ^,E n G Cn E n iP) + ^(V>, {E^ + E n ) 2 ^), 
2c 2 2 

where we used the inequality 

|(0i,^)|<^||0i|| 2 + p0 2 || 2 

in both terms of (27). Let 7 = min {7c ; } be the minimum 
of the gaps of the Hamiltonians ■ 

The first term on the rhs of (28) is less than 
2^(ip,hc n ip). Now, assuming the key estimate, 
||Ge„-En|| < e, we see that 



2 - ci - - ||£^|| 2 - C2 ||(£„-i + £„)V|| 2 

C17 



We sum over n using 
above and get 



2 _ 



= Eo<n<N \\EnW from 



2 - d - - - 2c 2 ) U\\ 2 - C2 ||^|| 2 < — (V, // [ll6] V). 
c 2 / C17 



Finally, we optimize the constants c\ , c 2 yielding ci = 
1 — e\/2, c 2 = e/\/2- This proves the gap inequality. □ 

In all the upcoming proofs on the various gaps we use 
the same definition of subsets C n ,A n of [1,6] and pro- 
jections G n ,E n . As usual y e [1,6] denotes the spot of 
the magnetic field. Let n;,n r be some non- negative in- 
tegers such that n r + ni > 1, and assume that n r > 0; 
the choice of ni, n r in general will depend on A and B. 
The idea behind the definition of C n is that we cover the 
chain [1,6] by adding points to an initially chosen inter- 
val Co = [y — ni,y + n r ] in an alternating manner. First 
we add a point to the right of Co, then to the left until 
we reach the point 1. Then we add points only to the 
right of C2( y -m-i) until we finish at b. More precisely, 
we define the sets C n in the following way: 

DEFINITION III.3. Let C = [y - n u y + n r ], where 
we may assume that y — nj — 1 < b — y — n r such 
that Co C [1,6]. The intervals for n > are then 
C\ = [y + n r ,y + n r + 1],C 2 = [y - n x - l,y - 
n{\, . . ■,C 2 (y-n l -i) = [1, 2],C 2 („_ ni )_i = [2y + n r - m - 
l,2y + n r - ni], . . . ,Cb-(m+n r )-i = [&-!,&]■ 

We start with the kink case. 

PROPOSITION III.4 (Kink sector, B\ + B\ > 0). 

Let tp be orthogonal to the ground state of the kink Hamil- 
tonian, H^~^(B), on a chain of length b. Then, there ex- 
ists a strictly positive function g + ~(B, A) and a function 
< e(-B, A) < l/\/2, which are both independent of b, 
such that the following gap inequality is satisfied 



(V, H+ M {B)^) > g+-(B, A)(l - V2e(B, A)) : 



(29) 



Proof. First we shift the ground state energy to be 0, and 
define the new Hamiltonian 

H [1M (B) = H+- ] (B)-^1. 

With the set-up from Definition (III. 3) we can write the 
Hamiltonian, H^ ^B), in the following form: 



b—(rn+n r ) — l 
i=0 



using 



H Co (B) for i = 
Kx+i for i > ' 



where Cj = [x,x+ 1] for some 1 < x < y — ni 01 y + n r < 
x < 6, and with h^~ +1 from (7). 

We can express the gap conditions as 



9 



1. H Aa (B) > gt; (B, A)(l - Go), where g+; (B, A)Js 
the gap for the finite chain Hamiltonian, H (B) . 

2. h+~ +1 = l-G[ xx+ i] for f < x < y-ni and y + n r < 
x < b. 

Let 7 = min {(7^~(B, A), 1} which is strictly positive. 
Finally, we need to verify the second condition in As- 
sumption III. 1 , and define 



\\G Cn n 



C n := sup 



So let tp satisfy 

G n tp = -0 and G n+ \ip = 0. 



(30) 



(31) 



First, let n — 2m, < m < y — ni — 1; the case n > 
2(y — m) is similar, and the case of odd 1 < n < 2(y — n;) 
will be considered later. 

Let tp be a ground state of A„, i.e. G n ip — tp such that 
Gn+iip — 0. Then with the definition from (14) 

y-\-n r -\-m 
i—y—rii—m 

where Xx( z ) i s perpendicular to Xx{z)- Let us make some 
definitions and call / := - ^b}-XbI ' X : = Xy+n r +m(z) = 
IT) + fq- nr - m \l), X 1 - ■= Xy+n r+m {z) = fq-"r-m-i 
|t) - ||>, and (1 + q 2 fl 2 \S) = q\M) - |tf>. Then, 



\GcM 



HICKCIx^)!! 2 
HxOx^ll 2 

q 2 1 + |/|2 g -2K+m+l) 

i + g 2 1 + 1/| V 2( "" +m) ' 



Let to = 0, then we choose n r such that rhs is less than 
1/2. The condition for G < 1/V2 is thus \f\q~ nr > 1. 
By monotonicity it is clear that the condition, C n =2m < 
1/V%, holds true for to > 0. 

Now we come to odd integers, n = 2m+l with < to < 
y — ni — 1. Let ip satisfy G n ip = ip, and G n +iip = 0, then 

V> is of the form $ = Xy_ n ,- m -i(*) ® 0;=^-™ X<(*)- 
We have thus 

IIGc^ll 2 



= 1 - 



1 1 + |/ , |2 (? 2(n 1 +m+l) 

1 + 9 2 l + |/|V (n,+m) ' 



Let m = 0, then we choose ni such that rhs is less than 
1/2. This is accomplished if 1 > \f\q ni . By monotonicity, 
C n= 2m+i < 1/ V2 for to > 0. Our condition for the choice 
of ni, n r is thus 



\B\\+Bo 



B 1 - iB 2 



< q 



-ni 



REMARK III. 5. It is clear that there always exist inte- 
gers ni and n r such that 1 > \f\q n ' and 1 < \f\q~ Hr 



\B\\ + B 3 



< 1, 



are satisfied. Now suppose that q < B iB 
then we may choose n/ = and n r = 1. In this case, 
g H (B, A) = g\ (B, A) is found explicitly in Appendix 
B, see (Bl). If B 3 = 0, then one needs to choose 
m = n r = 1 and diagonalize a three-site Hamiltonian 
which we will not do here. 

PROPOSITION III.6 (Droplet sector). B 2 + B 2 > 

0. Let ip be orthogonal to the ground state of the droplet 
Hamiltonian, H^^(B), on a chain of length b. Then, 

there exists a strictly positive function g ++ (B,A) and a 
positive function < e(B,A) < l/\/2, which are both 
independent of b, such that the following gap inequality 
is satisfied 

(<P,H++ b] (B)TP) > g++(B, A)(l - V2e(B,A)) 2 M 2 . 

(32) 

Proof. First, we need to shift the ground state energy, 
and define a new Hamiltonian 

H [1M (B) = H++(B) + j\\(Bi,B 2 ,B 3 A)\\ \A. 

Using the sets from Definition (III. 3) we have the decom- 
position 

6— (ri(+n r ) — 1 

H m(B) = £ h Ci , 

i=0 



with 



H Ca (B) for* = 

for some x: 1 < x < y — ni 
for some x: y + n r < x < b, 



li xx+l 

h — h 
"■xx+1 



depending on whether d is to left (right) of y. h^ x+1 and 
hxx+i are taken from equations (7), respectively (8). We 
have the following gap properties: 

1. H Ao (B) > g++(B, A)(l - Go), where g++(B) is 
the gap for the Hamiltonian, Hc {B). 



htx+i — 1 G[: 



c+i] for 1 < x < V - n i- 



3- ^xx + l — 1 



G^x+i] for y + n r <x<b. 



□ 



We are left with verifying the key estimate (25). So let 
7^ tp satisfy G n ip — tp such that G n +\ip = 0. If the in- 
terval Gj+i is to the left of y then we have the same situa- 
tion as in the previous proof with the condition 1 > \f\q n ' 
and the slightly modified / = ll(Bl ' B2 f i 3 ^ ) 2 ll+i?3 ~ A . 

If the interval Ci+\ is to the right of y, then we will 
arrive at the same condition for n r , namely 1 > \f\q Ur - 
This is true by symmetry but one can easily derive this 
in the very same way we did in the other case. □ 
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REMARK III. 7. The same remarks are in order here 
for the droplet Hamiltonian. So let us suppose that 



(B 1 ,B 2 ,B 3 -A)\\ + B 3 -A 



< 1, or cquivalcntly, A > B 3 , 
= and n r — 1. In this case, 



Bi-iS 2 

then wc choose n; 
g ++ (B,A) = g2 + (B,A) is explicitly calculated in Ap- 
pendix B, see (B2). 



PROPOSITION III.8 (Droplet sector). Let B x = 

B 2 = 0, and B < A. Let ip be orthogonal to the 
all spin-up ground state of the droplet Hamiltonian, 
(0, 0, B), on a chain of length b > 3, and let 

g 3 (B, A) be the gap for the three-site Hamiltonian from 
eq (B6). Then, 



i . H++{BW > 2g+(B, A) [ ± 



(33) 



Proof. Again, wc need to shift the ground state energy, 
and define a new Hamiltonian 

H [hb] (B)=H++ (0,0, B)-§1. 

As before, we use the same decomposition of H[iu(B) 
into local Hamiltonians, h^. The first gap condition of 
Co has to be changed into 

H Ao (B)>g+(B,A)(l-G ). 



We only need to compute 

C n := 



\\GcM 
sup I fn 

O^VeWA n+1 :BnV'=V' \W\\ 



So let us take a (non-zero) vector tp such that G n ip = ijj 
and G n+ \ip = 0. If C n+ i is to the left of y, then = | 1 
) ® |T •••?>, and G c „ = 1 — leXCI- Then' 

l|G C >|| 2 _ 1 



1 + a 2 ' 



which is less than 1/2, and (1 — \/~2e) 2 = 

2 (72 ~ ^qTq^) ■ 

By symmetry this is also the condition if C„+i is to 
the right of y. More precisely, V ;= |t"""T)®|l}> an d 
Gc B = 1 - |C) (CI with (1 + g 2 ) 1 /2|0 = gilt) - in). 

By choosing Cq — [y — l,y + 1], we have verified the 
statement. The three-site gap, g^(B,A) is calculated in 
Appendix B. □ 



below, that the product state ip(z — — ^L^^ Q y )i which 
is a ground state of Hq , is a also a ground state of 
V. Therefore, we have found a ground state of Hq (B). 
That it is the unique ground state follows by combin- 
ing two facts: 1) ip(z) is the unique kink state with this 
property, which we will show, and 2) the vectors ip(z), for 
arbitrary complex z, span the full ground state space of 
Hq 6 ' 22 and there is gap to the rest of the spectrum 23,28 . 
So, it only remains to prove that among all vectors ip(z), 
there is a unique one that is a ground state and that the 
corresponding value of z is as stated in the proposition. 

Since ip(z) is of product form, and V acts non-trivially 
only at site y, we are left to show that 



B- S y Xy(z) = -j\\B\\ Xv (z). 



(Al) 



Without loss of generality, we may assume that ||B|| = 1. 
Then, zq~ y = — gl^zff^ • Now, checking all 2j + 1 vector 
components in (Al) we obtain 

\pn{Bi + iB 2 )w n+1 { zq -yy- n - 1 + nB 3 w n ( zq -yy- n 

+\p n -i(B l - iB 2 )w n ^(zq-yy- n+1 
= JWn (zq-yy- n , 



with p n = + 1) — n — n 2 , and for \n\ < j. This 

leads to the following equation: 

ip„w„+i(l - Bl) + ip n _no n _i(l + B z f 

-nB 3 w n (l + B 3 ) 
= jw n (l + B 3 ). 

By a straightforward calculations one verifies that 
1 1 

-p„w n+ i + -p„-iw n -i = jw n , 

Pn-lWn-! = (j + n)w n , 

1 1 

- -^PnWn+l + -pn-lWn-l = nW n . 



This proves (Al). 



□ 



APPENDIX B: EXPLICIT DIAGONALIZATIONS 
OF SMALL-SITE SPIN 1/2 HAMILTONIANS 

1. H+-(B) 



APPENDIX A: PROOF OF PROPOSITION II.l 

Proof. First, it is clear that the bounded perturbation 
V = B ■ S y can shift the ground state energy of Hq 
by no more than its norm, j||B||. Wc claim, and show 



Here we diagonalize the two-site Hamiltonian, 
H^y (B), with magnetic field not parallel to z-axis at y = 
1. By the XX symmetry we may assume B 2 = 0. Since 
we already know two eigenvalues, namely, ±^\/Bf + B 3 , 
it is best to factor them out from the characteristic equa- 
tion. Another way is to diagonalize the Hamiltonian re- 
stricted to the orthogonal complement of the two known 
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eigenvectors. The Hamiltonian is of the form 
2H+~ ] (B 1 ,B 3 ) 



The gap above the ground state is therefore 



B 3 





B 1 








1 - A + B 3 


-A- 1 


Bi 


Bi 


-A- 1 


1 + A-B 3 








Si 





-B 3 



The characteristic polynomial, p, is equal to 

p(t) = (t 2 - Bl){t 2 - 2t + 2AB 3 - 2B 2 - B 2 ) 
+ 2B 2 (t - AB 3 ) + Bf. 

We divide this polynomial by t 2 — B\ — B\ (Note that 
we have multiplied the Hamiltonian by two) obtaining 

p(t) = (t 2 -2t-B 2 - B 2 + 2B 3 A)(t 2 - B 2 - B%). 

The two eigenvalues we are looking for are thus 

t± = 1 ± + Bj + B 2 - 2B 3 A. 
One can easily verify that 



B{ + B£ > 1 - yi + B( + B£ - 2B 3 A. 

Hence, the gap between the lowest eigenvalues of 
-ffjl 2] (Bi,B 2 , B 3 ) is equal to 



g+-(B) = \ \yj\ + \B\*- 2B 3 A+\\\B\\, (Bl) 
which is a positive function. 

2. H++(B),B! + Bl>0 

The diagonalization of the two-site droplet Hamilto- 
nian, H^(B), with the field at y = 1 is very similar to 
the two-site kink Hamiltonian. We have 

2H++(B 1 ,B 3 )-A1 

/ B 3 - A Si \ 

1 + B 3 -A- 1 B l 

B l -A- 1 l-B 3 

V B x A-B 3 J 

The characteristic polynomial, q, of the rhs is equal to 

q(t) = (t 2 - (B 3 - A) 2 )(t 2 - 2t - 2B\ - B\ + A 2 ) 
- 2Bf{t - A 2 ) + Bl 



t = ± \J B\ + (B 3 — A) 2 are two roots and we factor them 
out from q(t), and obtain 

q(t) = {t 2 -2t-B\-B\ + A 2 )(t 2 -Bj-(B 3 -A 2 )). 
The two new eigenvalues of H^(B) are 



'l + \\B\\ 2 - A 2 + A) . 



g+ + (B,A) = Ul-Jl + \\B\\ 2 -A 2 



+ \[b~1 + B 2 + (B 3 - A) 2 ^j . (B2) 

3. H++{Q,Q,B),B< A 

Since in this case there is only a magnetic field in the z- 
direction we can easily diagonalize the three-site droplet 
Hamiltonian, H = H^(0,0, B), with the field in the 
middle at y = 2. Then H commutes with the symmetry 
S : S(u ai Ua 2 Ucr 3 ) = u a3 u a2 u ai , where a l = ±. 
We choose the following eigenbasis of S: 



Vl = 


(1,0)®(1,0)®(1,0), 








V2 = 


(0,1)0(0,1)0(0,1), 








v 3 = 


^(1,0)0(1,0)0(0,1)- 


^(0,1) <s 


5(1,0)$ 


5(1,0), 


v 4 = 


^(1,0)0(0,1)0(0,1)- 


^(0,1)6 


5 (0,1) a 


5(1,0), 


V5 = 


^(1,0)0(1,0)0(0,1) + 


^(0,1)6 


5(1,0)$ 


5(1,0), 


v 6 = 


(1,0)0(0,1)0(1,0), 








V 7 = 


(0,1)0(1,0)0(0,1), 








v 8 = 


^(0,1)0(0,1)0(1,0) + 


^(1,0)6 


5(0,1)$ 


5(0,1). 



Vi, V2, v 3 , «4 are eigenvectors of H with eigenvalues ei = 
f , e 2 = A- § ,e 3 = \{A+l + B) and e 4 = \{A+l-B), 
respectively. What remains are two copies of the two- 
dimensional matrix (due to the symmetry S) 

N(B) 



1 / A + l + B -(V2A)- 1 

2 I -(V2A)- 1 A+l-B 



The matrix N is equal to H ++ {B) reduced to the 
span{w 5 , vq}, as well as to span{«7, v$}. The eigenval- 
ues are equal to 

e 5 = e 7 = \ {a + 1- ^\A- 2 + B 2 ^ , (B3) 

es = \[A+l + sJ\A- 2 + B^. (B4) 



Notice that for B < B = 3 \^if\ 
ei < e 5 = 67 < e 2 . 



(B5) 



This says that the lowest energies in the total S^-sectors 
are ordered (though not strictly) by their energy. The 
gap for B < A is equal to 



gt(B,A) 



+ A- y^A- 2 + B 2 - B I loi /;../; 

B 



(B6) 
for B < B < A 
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APPENDIX C: EXCITATION £-(B) 

Here we calculate the lowest eigenvalue of H(B) = 
+ + B(Sq — in the sector with one overturned spin. 
Since we want to avoid complications from finite chain 
boundary effects, we prefer to treat the infinite volume 
case with the magnetic field at 0, say. We first take B > 
0. 

Let tp = ^2 X a x S~ |t). Then, for tp being an eigenvector 
of H(B) with energy £, we have (| \S+H(B)\yb) = £a x , 
and thus we get the equations 

a x+1 = 2A(1 - £)a x - a x _i, for \x\ > 1 (CI) 
ai = 2A(l-f + B)a -a-i. (C2) 

It turns out that in addition to the pure absolutely con- 
tinuous spectrum of the discrete Laplacian (in the units 
here, it is the interval [1 — A -1 , 1 + A -1 ]) there are two 
(a highest and a lowest) eigenvalue generated by the per- 
turbation BSq. Let 

r± = A(l -£)± ^/A 2 (l-£) 2 - 1 

be the solutions to the characteristic polynomial. Then, 
all solutions of (CI) are of the form a x = ct\r x + a 2 r~ x 
for | a; | > 1. Notice that r + = rZ 1 ■ We now look for 
the solution a x = r' x ' with r = r_ which produces an 
eigenvector. With this choice, we have |r_| < 1, we insert 
this into (C2). Then we get 

AB = ^A 2 (l-£) 2 - 1, 



from which we conclude £±{B) = 1 ± y/B 2 + A~ 2 . From 
the gap at B = we know 23 that £(0) = 1 - A" 1 . 
Thus, the correct solution is £-{B) which, for the origi- 
nal Hamiltonian of interest, namely Hq + + BSq, has to 
be shifted back by B/2. 

Similarly, if B < 0, then we study H(B) = H+ + + 
B(Sq + \) which amounts to replacing B by — B in (Cl- 
C2). 

The lowest energy state of + + BS^ in the sector 
with one overturned spin is thus 

£-(B) = 1-Vb 2 + A- 2 + ^\B\. (C3) 
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